# Table S3: Descriptive statistics

# Load packages
library(psych)

# Import data
path <- "P:/2017-pathways/new/2-data"
setwd(path)
load("thetameta2.Rdata")

# describe data
myvars<-c("theta1", "theta2", "theta3", "theta4", "theta5", "theta6", 
          "theta7", "theta8", "theta9", "west", "female", "urban", "age",
          "treatment1", "treatment2", "treatment3", "he", "worried", 
          "lofoten", "oilwork2")
descdata<-thetameta2[myvars]
desc<-describe(descdata)
print(desc,digits=3)